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The type-I seesaw Lagrangian yields a non-generic set of active-sterile oscillation parameters — 
the neutrino mass eigenvalues and the physical elements of the full mixing matrix are entwined. 
For this reason one is able to, in principle, test the model by performing enough measurements 
which are sensitive to neutrino masses and lepton mixing. We point out that for light enough right- 
handed neutrino masses — less than 10 eV — next-generation short-baseline neutrino oscillation 
experiments may be able to unambiguously rule out (or "rule in") the low energy seesaw as the 
Lagrangian that describes neutrino masses. These types of searches are already under consideration 
in order to address the many anomalies from accelerator neutrino experiments (LSND, MiniBooNe), 
reactor neutrino experiments (the "reactor anomaly") and others. In order to test the low-energy 
seesaw, it is crucial to explore different oscillation channels, including Ve and ffj, disappearance and 
i^n i^T appearance. 

PACS numbers: 14.60.Pq, 14.60.St 

I. INTRODUCTION 

The simplest way of accommodating nonzero neutrino masses is to add to the standard model particle content at 
least two Weyl fermions Ni {i — 1, 2, . . .) which are not charged under the known SU{3)c x SU{2)l x U{1)y gauge 
group. The most general renormalizable Lagrangian is very well-known and easy to write down. Gauge singlet fermions 
— henceforth referred to as right-handed or sterile neutrinos — are only allowed to couple to the standard model via 
a Yukawa interaction with the Higgs and lepton SU (2) ^ doublets. On the other hand, right-handed neutrinos are 
allowed to have nonzero Majorana masses. In summary, the most general renormalizable Lagrangian is 

Mi 

- Com + Ckin - -Y^rN, ~ y'^'L^N.H -f H.c, (LI) 

where >Coid is the standard model Lagrangian in the absence of gauge singlet fermions, £kin are the kinetic energy 
terms for the Ni fields, a = e,/x,T are the flavor indices of the three lepton doublet fermion fields L, and H is the 
Higgs boson doublet field, y are the neutrino Yukawa couplings, while M are the right-handed neutrino Majorana 
mass parameters. Note that we picked a weak basis for the N states such that their Majorana mass matrix is diagonaL 

After electroweak symmetry breaking, Eq. (ILII) describes, in general, 3 + n (n is the number of right-handed 
neutrinos) Majorana fermions — neutrinos — most of which are massive if n > 2. It is very easy to accommodate 
the universally accepted neutrino oscillation data by judiciously choosing values for y and M. Assuming all M to be 
of the same order of magnitude, current data rule out 1 neV ^ M < 1 eV [H-lll . All other values are allowed! 

More information is required in order to test whether Eq. (jl.lj) is the correct way of understanding neutrino masses. 
Theoretical considerations allow one to rule out M > 10^^ GeV Q, while a simple interpretation of the gauge hierarchy 
problem leads one to favor M < 10'' GeV, assuming there are no other new states at or above the electroweak symmetry 
breaking scale @. Naturalness is not a good guide when it comes to picking a value for M — all values for M are 
technically natural in the sense that in the limit where all M vanish the global symmetries of Eq. (jl.l|) are augmented 
by U{1)l, global lepton number (see, for example, [L]).* 

Direct experimental probes of Ec^. p.ip are possible for M values below the TeV scale. Here we are mostly interested 
in M values in the few to several eV range. In this mass-range, right-handed neutrinos can be "seen" in neutrino 
oscillation data. Most interesting, in this mass-range, there is the possibility of unambiguously testing Eq. (jl.ip . The 
reason for this is that Eq. (11.11) implies relations among some of the mixing angles that can be measured in neutrino 
oscillation experiments, as we discuss in detail in Sec. |lTl In particular, given what is known about neutrino masses 
and lepton mixing today, simple lower bounds on different combinations of mixing parameters exist and, if these were 
to be violated by future experimental data, Eq. (11.11) would be ruled out for certain values of M. We discuss some of 
these quantitative bounds in Sec. IIIII 



* Lepton number, along with baryon number, is anomalous. Nonetheless, if n = 3, U{1)b—l is a proper non-anomalous global symmetry 
of Eq. ||LT] | in the limit M -i> 0. 
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Outside of M around a few to several eV, there are several interesting direct and indirect experimental bounds on 
Eq. (jl.ip which have been discussed in the literature. We briefly comment on these, and summarize our results in 
SeclWl 



II. FORMALISM AND RELATIONS AMONG ELEMENTS OF THE MIXING MATRIX 



After electroweak symmetry breaking, the (3 + n) x (3 + n) neutrino Majorana mass matrix from Eq. (jl.ip is 



03x3 ruD 
ni^ M 



(11.2) 



where O3X3 stands for a 3 x 3 zero matrix, M is the nxn matrix of right-handed neutrino Majorana mass parameters 
and Tod — yv is the 3 x n matrix of Dirac mass parameters, v is the vacuum expectation value of the neutral 
component of H. Eq. pi.2p is expressed in the so-called flavor basis: v^-^v^^Vt, i^si, ■ ■ ■ , i^s,,- This is related to the 
neutrino mass basis — z^i, 1^2, . • . , vs+n — via the (3 -I- n) x (3-1- n) neutrino unitary mixing matrix U: 

a = e, /X, T, Si, . . . , Sn, i ^ 1,2, . . . ,n + 3 such that 



, 03x3 run 
nl M 



U* 



I mi 



V 



TO2 



(II.3) 



where toi, m2, . . . , TO3+„ are the neutrino masses. It is instructive to express U in terms of 4 submatrices. 



(11.4) 



where Vmns is a 3 x 3 matrix that will be referred to as the "active" mixing matrix, l^g is an n x n matrix that will 
be referred to as the "sterile" mixing matrix, O, O' are 3 x n matrices and we will refer to as the "active-sterile" 
mixing matrix. In the absence of sterile neutrinos,^ Vmns is the standard neutrino mixing matrix whose elements have 
currently been measured with different degrees of precision ,7J. In practice, unless otherwise noted, we will assume 
that the elements of are parametrically smaller than those of Vmns and that Vmns is approximately unitary: 

VmnsK 



MNS 



L3x3- 



In the absence of interactions beyond those in Eq. (11.11) . physics is invariant under rotations among the sterile states, 



l3x3 03xn 
On X 3 X n 



u. 



(11.5) 



where 13x3 is the 3x3 identity matrix and Unxn is any nxn unitary matrix. This invariance allows one to significantly 
reduce the number of physical parameters in U and reveals that all physical phenomena are specified once the elements 
of ^MNS and O are known (these are invariant under Eq. (jll.5[) ). In summary, the "physical" parameters of U are Uai, 
where a = e, fj,,T, i = 1, 2, 3, . . . , 3 -I- n the active neutrino flavors. Usi, s = si, . . . , s„ cannot be directly measured. 
For more details on how to parameterize U in the case at hand see, for example, [8l-lll|. 

Eq. pi.3p further provides nontrivial relations among the elements of U and the neutrino mass eigenvalues. In 
detail. 



0i' 



(11.6) 



t In this case, of course, the neutrino masses would have to be generated through other means, or the neutrinos are Dirac fermions (all 
M = 0). 
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where a, /3 = e, /x, r, si, . . . , s„. We now proceed as follows: we will choose as inputs irii and the elements of Vmns 
and try to predict the elements of 9. All the extra information comes from the vanishing "active-active" elements of 
m„ — 03x3: 

J2t UaiUurrii = 0, 

SLl UaiUb^m^ + J27=l Ua,i+3Ub.i+3'mi+3 = 0, 
J27=l'^a.iQbim.i+3 = - I]Ll(^MNs)ai(VMNs)6jTOi, (11.7) 

for all a,b — e, /i, r. Note that we choose all nii to be real in order to re-express Eq. (|IL6p in terms of the elements of 
U, as opposed to U*. 

A few simplified scenarios illustrate the spirit of the results we will discuss below. For example, if there was only 
one active Va and one right-handed neutrino Vs, Eq. (jll.7|) would read 

C/,>4 = -U^i^m, {11.8) 

where, in order to facilitate future comparison with more complicated scenarios, we label our masses mi and 7714. In 
this case U can be chosen real and, if we further redefine the neutrino mass eigenstates in such a way that mi — —mi, 
it is easy to see that the only parameter^ necessary to parameterize U is uniquely determined by the neutrino mass 
eigenvalues: tan^ = mi/m^. This scenario could be tested if one were to independently measure the two neutrino 
masses and the one leptonic mixing angle. 

In the case of one active and two sterile states, Eq. pi. 71) reads (after redefining toi — > —mi) 

U^^mi + Ul^m^ = f/^imi . (II.9) 

In this case, it is possible to obtain lower bounds for the active-sterile mixing angles as a function of Uai and the 
masses. A formal solution to the constraint above is 

C/a^4 = C/fl— C0S\, (11.10) 

m4 

^f5 = C/,^i— sin^C (11.11) 
ms 

where C € C. One can choose Uai real and easily show 

\Uai\^mi + \Ua5\hn5 > \Uai\^mi. (11.12) 

The inequality is saturated for 5(C) — 0. If either entry of the active-sterile part of the mixing matrix vanishes, 
the magnitude of the other one is guaranteed to be equal to C^aiy^^^^- The careful reader will note that unitarity 

imposes other constraints, including, |?7qiP + |?7a4p + |C^a5p = 1. As mentioned in passing earlier, we are working 
under the assumption that m4^5^... ^ mi^2,3 and the elements of Vmns are parametrically larger than those of Q. In 
this case, this assumption translates into Uai = 1 + 0{mi/mm^ 5)^. 

In the case of two active states and i^b and one sterile state, Eq. (jll.7p reads (after redefining m4 — > — m4) 

U!^m^ = U!^m2, (11.13) 
[7,^4 = (11.14) 
Ua4Ub4m4 = Ua2Ub2m2, (11.15) 

where we took into account that, in this case, one of the neutrino masses (here mi) vanishes. The three equations 
above translate into 

(11.16) 
(11.17) 

Here, all elements of U can be chosen real. Requiring m4 ^ m2 we can approximate 11^2 = 1 — Ub2 = sin^ iJ, and it 
is easy to see that the active-sterile entries of the mixing matrix are well-define functions of the active mixing angle 





-^"2m4' 








— '-'b2 




m4 



t We choose the standard parameterization for the 2x2 mixing matrix, Uai = Us4, = cos 6, Ua2 = 



—Uai = sin 9. 
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•0 and the neutrino masses, 7712,7714. By measuring "active-active" oscillations governed by the small Aml2 = 77i| 
mass-squared difference one can in principle measure i? and rn^ and hence express the oscillation probabilities at short 
distances as a function of only one new parameter (1714). For example, for an ultrarelativistic beam of neutrinos with 
energy E and short baselines L {ra\L/E <^ 1), 

P,,^l-4sin2z?(l^)sin2(^), (11.18) 

P,t-1- 4cos2 ^ (h^) sin^ (^) , (11.19) 

P.,~sin2 279(ll^)%in^(!#). (11.20) 

For more complicated cases, including the physically relevant case of three active neutrino flavors, it is convenient 
to formally solve Eq. (jll.7p . which can first be re-expressed as a matrix equation. After redefining the sign of the light 
neutrino masses (mi — >■ —mi, etc), Eq. (jll.Tp reads 

( 6 ^TTlhcavy ) ( 6 V^hoavy ) ^ = ( VmNS V^ligh* ) (^MNS y/might)'^ , (11.21) 

where TTiheavy = diag(r774, . . . , 7773+„), miight = diag(mi, 7712, m^). The equation above allows one to write 

6V"^hcavy = VmNS V™light-R, (11.22) 

where i? is a complex matrix. This very useful solution was flrst presented in [l2j . In the case 77 = 3, i? is a complex 
orthogonal matrix, RR^ ~ R^R ~ 1. Here, we will also be interested in the case 77 = 2 when R "contains" a complex 
orthogonal matrix,^ 

/ \ / cosC sinC\ 

R3X2 — cos( sinC or — sinC^ cos^ . (11.23) 
y -sinC cosC J \ J 

We choose the vanishing neutrino mass eigenvalue to be mi = or 7713 = 0, respectively, and ^ is a complex number. 
We draw attention to the decisive role played by the zero eigenvalue in rendering Eq. pi. 221) a solution to Eq. pi.2ip . 

Before moving on to the physically relevant cases, we quickly discuss the case of two active states and two sterile 
states. Eq. pl.22p reads 

/ UaAy/ml Uab\/m^\ _ ( COS ■dy^mie^'^ sin 79.^7712 \ / cos^ sin^' 
\^ ?7fc4V7774 Ub5y/rn5 J y — sin 79^77716*'^ cos79y^m2 J \^ --sin^ cos( 

where C S C and we assumed 7714.5 3> 70,1^2, and the active-active mixing matrix is approximately unitary and 
parameterized by one mixing angle 7? and one "Majorana" phase (p. Overall, the 4x4 mixing matrix is parameterized 
in terms of four real parameters, 7^, 1^9, and the real and imaginary parts of C.. 

In order to separate the impact of the relations among the different elements of U on the possible values of active- 
sterile "mixing angles," it is convenient to work with the dimensionless objects Xaj = Uajy/mj/y/m^ [j = 4,5) and 
c = ^mij ^m^e^"^ , where m-i is deflned as the largest of m\^2- Hence 

/ Xa4 Xah \ _ ( c cos 7? cos ^ — siu 7? siu ( c cos 7? sin -|- sin 7? cos C 
y Xb4 j \^ — csiuT^cosC — C0S7?sinC — c sin 7? sin ^ -I- cos 7? cos C 

In general, one can always pick C such that one of the Xaj vanishes. When that happens, however, none of the 
other three vanish.^ For example, Xas, = implies tan^' — cj tan 7?. Under these circumstances, all other elements of 
Xaj are uniquely determined functions of the three (real) active neutrino parameters, 79,c. 



(11.24) 



(H.25) 



§ In the "1+2" case discussed above, Rix'2 = (cosi^ sin(^) and RR?" = 1, while in the "2+1" case i?,2xi = (0 1)"^. In the latter case, 
iJi?^ = diag(0, 1). 

^ We are ignoring special values of, or special relations between, •& and c. These are, after all, independent, a priori uncorrelated, 
parameters. 
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III. QUANTITATIVE UPPER BOUNDS 



Here we quantitatively discuss the expectations for the values of the elements of in the physically interesting case 
of three active neutrinos and n = 2 or n = 3. The n = 2 case is both simpler to analyze and "more minimal," so we 
devote most of our presentation to it. We first quickly summarize some of the inputs we will use throughout. 



A. Active Oscillation Parameters 



Current experimental data constrain all the elements of Vmns and two mass-squared differences, 



Throughout we assume sin 6*12 = 0.3, sin^6': 



23 



0.5. 



Am^2 and 



7.6 X 10-^ eV^ and jAmfg] = 2.4 



agreement with recent analysis of all neutrino data [ij; [3 ■ Unless otherwise noted 
In summary, we assume 



_2 

X 10^ 



' eV^ in 
we will assume sin'^ = 0.01. 



0.83 


0.55 


0.1e"5 \ 









-0.39 - 0.06e-"' 


0.59 - 0.04e^*'' 


0.70 


U 


e# 




0.39 - O.OGe-"' 


-0.59 - 0.04e-*'' 


0.70 / 










(III.l) 



The CP-odd phase 5 is currently only very poorly constrained and will be taken as a free parameter when relevant. 
Note that, since ^13 is known to be small, the impact of the uncertainty in i5 on the magnitude of the elements of 
Vmns is relatively small. The so-called Majorana phases tp, (p will be taken as free parameters when relevant. When 
one of the mass eigenvalues vanishes, only one Majorana phase is physical.^ 

In the case of a so-called normal neutrino mass hierarchy, we express all neutrino masses as functions of mi: 



771, = 777f + 



At 



'12; 



7775 = 7777 



Note that all masses are defined as real and positive. The square-root of the 



matrix 777iight is 



V 



/lEi 

rn 









1/ 



7771 



= 0.22 




(III.2) 



In the case of a so-called inverted neutrino mass hierarchy, we express all neutrino masses as functions of 7773: 



Ar 



n3 



and 



Ar 



+ Ar 






^12- 



We remind readers that A777 i3 is, in this case, negative. The 



7773 



^Tm^=0.22 




(III.3) 



square-root of the matrix 7T7iight is 

/ 

V 

Similar to what was done in the end of the last section, we will be mostly interested in computing bounds that are 
independent from the values of the heavy masses. With that in mind, we define the matrix X = G^mhcavy/ ^7773 
(normal hierarchy) 01 X = 0.y777hcavy/ -\/"^2 (inverted hierarchy). It is trivial to convert between a bound on Xai and 
Uai- For example, for a normal hierarchy, Ue5 = X^z y^ms/m^^ and, for an inverted hierarchy, Ut4 — Xr4 -yAn^/rol . 

With neutrino oscillation experiments one can study, in principle, the disappearance probability of active neutri- 
nos, Paa, or the conversion probability of one active flavor into another, Pat {a,h = e, /i,t). These, in turn, are 
proportional to mass-squared differences and products of elements of the mixing matrix. Here, we will be interested 
in next-generation short baseline experiments. Short-baseline, here, means values of L such that |A777^3|L/£' ^ 1, 
so "atmospheric" induced oscillations (and, hence, "solar" induced oscillations as well) can be safely ignored. Under 
these circumstances one is only sensitive to the large mass-squared differences. In the limit r77hcavy ^ "T-iight, at short 
baselines, we will assume one is able to measure all different oscillation frequencies, proportional to 777^, j — 4,5,6, 
and the following combinations of mixing elements: |?7ajP, j =4,5,6 (from disappearance) UajU^p j = 4,5,6 (from 
appearance).^ Henceforth, unless otherwise noted, we will concentrate on placing bounds on l^ajp and for 



* We pick the "standard" parameterization for the active-active mixing matrix and the mass-squared differences, along with the "standard" 

definition of the neutrino mass eigenvalues. For more details see [H and the many references therein, 
t ff mi = 0, can be chosen as the physical Majorana phase, fn the case 7713 = 0, one can choose </> — ■(/> as the physical Majorana phase, 
•t it is convenient to label the elements of Xai such that a = e, ^, r and jf = 4, 5, . . . 3 + n. We do that henceforth. 
§ Needless to say, experimentally, this is a formidable task at best! 
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different a,b = e, ^, r. 

Before proceeding, we repeat our conventions and the approximations that go into our resuhs. We assume there 
are three active and n = 2 or 3 sterile neutrinos. After electroweak symmetry breaking, the neutrinos "spht" 
into three hght, mostly active neutrinos with masses 7711,7712,7713 and 77 heavy, mostly sterile neutrinos with masses 
777,4, ^775, . . . , ms^n- We also assume that Vmns is approximately unitary and that its elements are known, except for the 
CP-odd phases ("Dirac" and "Majorana" ) . As far as the "scale-free" active-sterile mixing matrix X is concerned, our 
approximations are safe if we restrict all Xaj ^ 1. This way, all elements are constrained to be less than or of order 
'T7iight/'77hcavy, which is Small as long as we stick to 777heavy values above 1 eV and stay away from quasi-degenerate 
values for mi, 7n.2, 7773. 



In the case of only two right-handed neutrinos, the lightest neutrino mass is zero and the R matrix is given by 
Eq. (III.23p . For a normal neutrino mass hierarchy. 



0.236'-^ O-le'-^ 
-(0.25-h0.02e-^*)e*'^ 0.70 



X™i=l (0.25 - 0.02e-V^ 0.70 I ( 1^ ) ■ (III.4) 



— sin C cos C 



On the other hand, for an inverted mass hierarchy. 



Xi„verted= I " (0.39 + 0.06e-'*)e*'^ 0.59- 0.04e-'^ I ( _'=°f^^ ). (III.5) 



0.83e*'^ 0.55 
y^.^J + Ome~'^)e'''' 0.59-0.0,^ M _ • r r 
(0.39 - 0.06e-*'')e''^ -0.59 - 0.04e-''5 / ^ 

Before discussing bounds on the different elements of X, we discuss two interesting and potentially relevant examples. 
For C = 3/47r + i, 6 = 6/57r, (/) ~ 7r/2 and a normal mass hierarchy, 

/ 0.41e-° *^6i o.45ei "3' \ 

^normal = OMe^'^^' O.eie-^'^^. (jjl.6) 

\ 1.27e2-44' 1.26e-2 «» / 

On the other hand, for ^ = 2/37r + 0.37, 6 = 0, tp — 7r/2, and an inverted mass hierarchy, 

/ 0.44e-2-24i o.62ei *^3i \ 

^inverted = 0.69e2-66» O.GGe-^^^ . (III.7) 

\0.71e-0-39' 0.60e"-89' / 

Both of these are consistent with the best "3-1-2" fits to the so-called short-baseline anomalies from LSND [l5| . 
MiniBooNE (T6l - [l8| and reactor data [l^. From [23|, for 7714 = 0.5 eV^ and 7^75 = 0.9 eV^, the best fit point occurs 
when |A:e4| = 0.49, \Xe5\ = 0.62, 1X^4} = 0.65, jAT^sl = 0.66 and arg(A:*4A:^4A:*5A:^5) = 5.0. These are qualitatively 
in agreement with the matrices above. ^ In the case of a normal mass hierarchy, a good fit can only be obtained if 
has an order one imaginary part. For real |A'e4^5| are always too small. In the case of an inverted hierarchy, one 
can qualitatively "fit" the short-baseline anomalies even for real C- We note that a qualitatively similar analysis of 
the short-baseline anomalies was presented in [2l|. 

Moving away from trying to explain the short-baseline anomalies, it is easy to see that, in general, for either mass 
hierarchy, at most one Xaj can vanish.** The situation here is similar to the 2-1-2 case discussed in the previous 
section. This means that if some \Xaj \ are constrained to be very small, the scenarios under consideration will predict 
that other \Xaj \ are larger than some lower bound. It also implies that if one can constrain two different \Xaj\ to be 
smaller than some amount, the model can be conclusively ruled out. 

We begin by discussing the most conservative scenario: 777.4 ^ ^5- In this case, all practical information concerning 
testing the model will be provided by Xa4, a = e, /j,, t. When 777.4 ^ '775, all Ua5 are parametrically much smaller than 



^ The careful reader will note that the relevant relative phases in the examples are on the small side. 

** Exceptions exist, of course, but these are entirely driven by the specific values of entries of Vmns- For example, for Ue-i = and 
Ut3 = U/iS (vanishing 613, maximal "atmospheric" mixing), it is easy to see that, for an inverted mass hierarchy, |X^4_5| = |Xr4,5| such 
that if one picks a value of ( that leads to a vanishing X^i, X^-n would also automatically vanish. 



0.05 0.1 0.15 0.2 0.2 0.4 0.6 




0.05 0.1 0.15 0.2 



|Xe4l 

FIG. 1: Boundary of the lowest allowed values of mixing parameters in different \Xa4\ x \Xb4 \ planes, a,b = e, /i, r, in the case 
of a normal neutrino mass hierarchy and n = 2. Different X are defined in Eq. (|TlL4)) . The regions located to the left/bottom 
of the curves are not accessible. 

Uai (by a factor mn/m^) and the only way to unambiguously test this model is to look for neutrino oscillations 
governed by the m4-induced oscillation length. Figures[I]and[2]depict, for a normal and inverted hierarchy, respectively, 
the lowest accessible values of pairs of |Xe4|, \X^4\ and |Xt-4| once one allows for all possible values^^ of C, <5, 4> or i^, 
respectively, as defined in Eqs. (jIII.4IIII.5[) . The figures allow one to quickly confirm some of the statements made 
above. At most one of |Xa4| vanishes for any point in the parameter space and there are "forbidden regions" (here we 
highlight the ones in the lower-left-hand side of the curves) in the space of the |v'Cq4|. If neutrino oscillation data were 
to limit these mixing parameters to within these forbidden regions, we would be able to rule out the model. It should 
be clear that by 'model' we refer not only to a low-energy seesaw Lagrangian with two right-handed neutrinos, but 
also to the mass of the lightest of the mostly sterile states, in this case TO4. Next, we discuss two concrete examples 
(i.e., constraints for fixed values of 7714). 

When m4 = 1 eV, current electron neutrino disappearance data constrain |Xe4| < 0.45 at the 90% confidence 
level (C.L.) [l^, while muon neutrino disappearance data constrain X^^ < 0.64 at the 90% C.L. [2^ (see also [ll])- 
Furthermore, searches for — )■ reveal that |Xe4||X^4| < 0.46 at the 90% C.L. [l^. Figs. [T] and [2] reveal that the 
n = 2 scenarios are consistent with these results for either neutrino mass hierarchy. 

In the case of a normal neutrino mass hierarchy, the current data allow either \Xf,4\ or |Xp4| to be vanishingly small, 
in which case the other mixing element (|^;i4| or |^e4|) is close to current experimental constraints. In the case of an 
inverted neutrino mass hierarchy, the situation is much more constrained. The constraints spelled out above translate 
into, at the 90% C.L., \X(.i\ > 0.34, \X^4\ > 0.60, and \Xr4\ > 0.61. These, in turn, allow one to make very specific 
predictions concerning short-baseline oscillation experiments sensitive to a new mass-squared difference aroud 1 eV^ . 
These include the observation of . . . 

• disappearance with an associated effective mixing angle sin^ 2-&ee > 0.02. An interesting new proposal to 



We numerically restrict the magnitude of the imaginary part of ^ during our exploration of the parameter space in order to satisfy the 
constraint \Xa-i\ < 1. 
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closely expose the Daya Bay detectors to a strong /3-eniitting source would be sensitive to sin^ 2dee ^ 0.04 [25l |: 

• disappearance with an associated effective mixing angle sin^ 2^?^^ > 0.07, very close to the most recent 
MINOS lower bound; 

• -;-> i^e transitions with an associated effective mixing angle sin^ z?e^ > 0.0004; 

• i^fj, i^T transitions with an associated effective mixing angle sin^ iJ^^. > 0.001. A ^ v^- appearance search 
sensitive to probabilities larger than 0.1% for a mass-squared difference of 1 eV^ would definitively rule out 
rui = 1 eV if the neutrino mass hierarchy is inverted. 

When 7714 = 10 eV, current electron neutrino disappearance data constrain \Xei\ < 2.85 at the 90% C.L. [26j . 
while muon neutrino disappearance data constrain |X^4| < 1.01 at the 90% C.L. |27[. Furthermore, searches for 
fij, — » z^e reveal that |Xe4||X^t4| < 2.9 [2^ and those for — > reveal that 1X^,411X1-41 < 1.8, all at the 90% C.L. 
[29| . Figs. [Hand reveal that the n = 2 scenarios are consistent with these results for either neutrino mass hierarchy. 
It is not surprising that the bounds on Xai are weaker than those for — 1 eV given that |C/a4p oc \Xa4\'^ /ni4. For 
m4 = 10 eV and all |Xa4| = 0.5, the related effective mixing angles are sin^ 2-daa = 0.005 (appearance experiments) 
and sin^ 2'dab = 6xl0~^. In order to unambiguously rule out 7714 = 10 eV and tt, = 2, one needs to either probe Ve 
or disappearance at the sub-percent level or look for appearance processes at the several parts-per-million level. 
Both appear very challenging, but the possibility of very sensitive — >■ v^- searches is currently under investigation 
(see, for example, f30|). 

In the less conservative case where both 7714 and are "small" (say, below tens of eV) both new oscillation 
frequencies are accessible to searches for neutrino oscillations at short baselines and the n = 2 low-energy seesaw 
scenario is much more constrained. Figure |3] depicts, for a normal and inverted hierarchy, the lowest accessible values 
of the magnitude of the active-sterile entries of the neutrino mixing matrix in the |Xa4| x {Xa^l plane, as defined in 
Eqs. (jIII.4iIII.5p . In our scan of the parameter space, we restricted the imaginary part of C to G(C) G !]• As 
advertised earlier, only at most one of the Xaj can vanish for j ~ 4,5 and a = e, fi,T. 
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FIG. 3; Boundary of the lowest allowed values of mixing parameters in different |Xa4| x \Xa5\ planes, a — e, /i, r, in the case 
of a normal (top) or inverted (bottom) neutrino mass hierarchy, and n — 2. Different X are defined in Eqs. (|III.4IIII.5|) . The 
regions located to the left/bottom of the curves are not accessible. 



The curves in Fig. [3] are easy to understand qualitatively. Eq. (|II.7|) , for a = b and n = 2 implies 

I^a4p + |^a5|' >a, (III.8) 

where Ca are positive constants that depend only on active neutrino parameters. We will further explore this inequality 
in the next subsection. Hence one expects the boundaries in Fig. [3] to resemble quarter-circles. This is not the case 
for a = and an inverted neutrino mass hierarchy. The reason is that in order to saturate the bound in Eq. (jIII.8|) . 
the magnitudes of some of the other Xaj grow too large, violating the assumptions we made in the beginning on this 
section. 

Fig.|3]also reveals that, in the limit where the magnitudes of the Xa4S are small, those of the Xa^^s are "large." This 
implies that the bounds discussed in the ^ tos limit above are more stringent if driven oscillations are also 
accessible to neutrino oscillation experiments. It is, unfortunately, cumbersome and less than illuminating to discuss 
bounds on the low-energy seesaw with two light right-handed neutrinos. The reasons are simple to understand. For 
example, one needs to discuss constraints as a function of both m4 and m^. Furthermore, oscillation lengths depend 
not only on Amf^ and Amfg, i — 1,2,3 but also on Am^g, while the oscillation amplitudes will depend on several 
different combinations of the J7a4,5, including relative phases. This all implies that, even though there are upper 
bounds on the magnitudes and relative phases of the elements of Ua4,5, they depend non-trivially on 7714 and ms. 

Here, instead, we will discuss the limiting case where 7714 and 7715 are very similar. In this case, all short-baseline 
oscillation probabilities will depend, approximately, on effective mixing angles, given by 



sm' 



2 27?„„=4(|t/a4p + |(7a5p), (HI.9) 



Sin^ 2dah = 4 \UaiU^i + Ua5U;^\^ . (IH.IO) 

where a,b = e, fj,,T and a ^ b. It is easy to see from Fig. |3] that none of the -daa can vanish, regardless of the mass 
hierarchy. Numerically, defining = ni^ = M, expectations for disappearance experiments are 

1 eV 1 eV 
sin^ 2i9ee > 0.009 x or > 0.08 x , (IH.ll) 

1 eV 1 eV 

sin^ 2d^^ > 0.08 x or > 0.02 x — — , (IH.12) 

, 1 eV 1 eV , , 

sin^ 2-drr > 0.08 x or > 0.02 x — — , III.13) 

for a normal or inverted neutrino mass hierarchy. Keep in mind that, as in the ^ case discussed above, these 
bounds cannot all be saturated at the same time. Hence, for M values less than 10 eV, one should observe i'^ (i^e) 
disappearance at more than the one percent level in the case of a normal (inverted) neutrino mass hierarchy. 
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In the case of appearance experiments, the foUowing upper bound apphes: 



sm^2'df,r > 1.7 X 10" 



1 cV 
M 



or > 3 X 10"^ X 



leV 

M 



(III.14) 



for a normal or inverted neutrino mass hierarchy. No nontrivial bound exists for sin^ ^-defi and sin^ ider . AU lower 
bounds are consistent with the current short baseline data for any M > 1 eV. Expectations for ly^ — > Vr appearance 
would be challenged by next-generation searches for short baseline Vt appearance from conventional neutrino super- 
beams for M values less than several eVs, for either mass hierarchy. As before, the different inequalities cannot all 
be saturated at once so the combination of results from different searches is more restrictive than the bounds of each 
individual search. 



In the case of three right-handed neutrinos, X is a 3 x 3 matrix. Compared to the n ~ 2 case, the number of physical 
parameters we need to consider is much larger. As far as the "active" parameters are concerned, the lightest neutrino 
mass is no longer zero (but constrained to be less than 0.3 eV or so [7]) and there are two physical Majorana phases, 
tp and (p, regardless of the mass hierarchy. The R matrix (see Eq. (jll.221) ') is now a complex 3x3 orthogonal matrix 
parameterized by 3 complex angles (here CtV-iO^ six real parameters. The analysis of the allowed values of X is 
much more cumbersome, and the results one can obtain are much less constrained. For this reason, our discussion 
will be briefer and less quantitative. Nonetheless, we will offer some qualitative statements that will prove useful to 
understanding the potential impact of next-generation experiments and how they might test the n = 3 scenario in 
the future. 

In the case <^ m5,me, neutrino oscillation experiments are only sensitive to oscillations involving the z/4 state. 
The relevant entries in X can be written as 



Xe4 = (VA/[Ns)elA/ — COS C COS ^ -|- (VMNs)e2W — Siu C COS ^ + (VMNs)e3 siu^, 

V m3 V "^3 

X^4 = (Vmns)miJ — - cosCcos^ -I- (yMNs)Ai2W — sin ( COS ^ -|- (Vmns)m3 sin^, 

V 7713 y 7773 

Xr4 = (^MNs)rlW — COSCCOS^ -I- (14lNs)r2 W — siu C COS f -|- (14lNs)r3 siu ^, 

V "73 V "^3 

in the case of a normal neutrino mass hierarchy, and 

Xe4 = (VMNs)elJ — COS ( COS ^ -|- (VMNs)e2 siu ( COS ^ -|- {VMNs)e3\ — siu^, 

V 7772 y m2 



X 



(VMNs)Atll 



■ cosCcos^ + (Vmns)^2 sin C cos ^ -I- (Vmns)a 



■sin^. 



(III.15) 
(III.16) 
(III. 17) 

(III.18) 
(III.19) 



Xr4 = (Vi 



MNSJtIi 



■ COSCCOS^ + (VMNs)r2 siu C COS ^ [V^ 



MNSJt3i 



■ sin^, 



(III.20) 



in the case of an inverted one. In either case, it is clear that one can pick values for the complex angles C and ^ 
such that two of the three Xa4 are vanishingly small. This renders this scenario impossible to definitively rule out in 
practice for any value of 7714 (as long as T4ins is unitary "enough"). If C, and ^ are such that Xf.4 and X^4 vanish, the 
only way to constrain this model would be to search for Vj- disappearance. Needless to say, experimentally, we are 
very far away from being able to do that with any precision. 

It is also possible to choose parameters such that all Xa4 are very small. In the case of a normal hierarchy, this can 
be achieved if 7711 <C 7772 and both sin^ and sin^ are very small in magnitude. In the case of an inverted hierarchy, 
the same effect can be obtained if 7773 ^ 7772 and | cos^| <C 1. Under these circumstances, qualitatively, the lightest 
mostly sterile neutrino "couples" predominantly with the lightest mostly active neutrino {vi or v^, depending on the 
mass hierarchy) and effectively decouples from the observable world in the limit where the lightest mass — currently 
unconstrained by experiment — is very small. In the limit that the lightest mass vanishes and this state becomes 
completely invisible,"''"'' we are back to the 77 = 2 case where plays the role of and vq that of v^. Very small 



tt In the case of an inverted mass hierarchy, this would happen when cos § = and ma = 0. 
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FIG. 4: Lower bounds for sm^ 2^ee x (M/1 eV) (black, solid line), sin^ 2i9p^ x (M/1 eV) (red, dashed line), and sin^ 2i9rT x 
(Af/1 eV) (blue, dotted line) as a function of the smallest neutrino mass mughteat for a normal (TOP) and inverted mass 
hierarchy (BOTTOM), for 71 = 3 and m4 = ms = me = M. 



values for all Xai are required if one wants the lightest mostly sterile state to play the role of warm dark matter, as 
discussed in j311] . 

In the case TO4 ^ ^ mg, the situation is much more constrained. As in the n = 2 case, we will restrict our 
attention to effective disappearance and appearance mixing parameters, 

Sin2 2^aa = 4 (|C/a4p + |C/a5p + I C/a6 P) , (111.21) 
Sin2 2^ab = 4 paiU^i + UaroU^s + UasU^f . (III.22) 

We performed a scan of the unconstrained part of the parameters space. In the n — 3 case, there are ten parameters: 
three complex mixing angles in R (C, '7,0' Dirac and two Majorana phases in Vmns f-i^d the smallest 

neutrino mass (toi or ma for a normal or inverted hierarchy) while keeping the other parameters fixed at their 
predetermined values, as described in Eq. (jlll.ll) and the text that surrounds it. Figure |4] depicts lower bounds 
for sin^ 2'daa x (M/1 eV) as a function of the smallest neutrino mass mughtest^^ for both normal and inverted mass 
hierarchies. Similar upper bounds for sin^ 2'dab prove to be less useful. The reason is that cancellations among 
the three different contributions to sin^ 2'dab allow for very small lower bounds that are, for all practical purposes, 
vanishingly small. 

Figure HI reveals that, for most values of ^lightest , the lower bounds here are similar to the ones obtained in the n — 2 
case when = m^. In fact, for mughtost ^ 10~^ eV, values much smaller than those of the other neutrino masses, 
the lower bounds are virtually the same. For intermediate values of rTiughtost, "cancellations" can occur. These, it 
turns out, are somewhat familiar to those well- versed in the neutrino physics literature. 

The lower bounds depicted in Figure |3] are very simple to understand. When = = mg = M, Eq. (jll.Tp . for 



"^lightest = ™i in the case of a normal hierarchy, mughteat = '"^S in the case of an inverted one. 
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a = b, translates into 



3 

y~^(VMNs)ai(^MNs) 



ai 



mi = m, 



aa- 



(III.23) 



Keep in mind that the different terms above are, in 



g' 



;eneral, complex. This, in turn, implies 



4(|[/a4p + |C/a5|' + |C/|^6) 



= Sin^ 2daa > 4 




(III.24) 



One way to saturate the bound is to choose ^, ry, C. such that two of the J7a4.5,6 vanish. Some care is required when 
niaa vanishes exactly (which happens for a range of values of rniightost for either hierarchy, see Figure |4]) , but it is 
safe to state that, under this circumstance, the lower bounds for some sin^ 2daa are vanishingly small for all practical 
purposes. The same argument also holds in the n — 2 case, as briefly discussed in the previous subsection. There, 
none of the \maa\ values are allowed to vanish. 

niaa, curiously enough, but not surprisingly, is the magnitude of the aa-element of what can be refereed to as 
the "active-active neutrino mass matrix," ruab- "m-ab is the effective neutrino mass matrix in the limit of very heavy 
sterile neutrino masses, and is also the active Majorana neutrino mass matrix that is generated in the so-called 
type-II seesaw (in a nutshell, the Standard Model Lagrangian augmented by an SU(2) triplet scalar that acquires a 
vacuum expectation value). The attentive reader will have noticed that the lower bound curve for sin^ 29ee is exactly 
proportional to the lower bounds on the magnitude of (often referred to as mfl« in the literature), used to estimate 
the sensitivity of neutrinoless double-beta decay experiments (see, for example, [7|). The equivalent observable with 
muon neutrinos, the magnitude of m^^, is less familiar but has been studied in the literature along with the magnitude 
of the other elements of niab [111- In the Appendix, we plot the allowed values of the magnitude of m^^, for both 
mass hierarchies, as a function of mughtcstj for the values of the mixing parameters chosen here. It is curious that 
vanishes only for an inverted mass hierarchy, and only in the limit when ~ m2 (to be more precise, > TO2/3). 
This is very easy to understand from the definition of m^^ and the known values of the active neutrino parameters 



Finally, in the n = 3 case, one may also consider the case m4 < <C rriQ. Qualitatively, we anticipate that 
the situation will resemble the n = 2 case, especially in the limit where we treat and V5 as one "effective" state. 
Scenarios of this type may be phenomenologically interesting. For example, V4 and may provide an explanation to 
the reactor anomaly, while ve may play the role of warm dark matter, as discussed, for example, in 



We have explored the seesaw Lagrangian (Eq. p.ip in the limit Mi ^ y^^v), concentrating on the values of the 
active-sterile neutrino mixing parameters. The key point is that even though the neutrino mixing matrix U is 
(3 -|- n) X (3 -f n) and unitary, it is not a generic unitary matrix. Several of its parameters are not observable and can 
be "rotated away," while the remaining parameters are related to one another and to the neutrino mass eigenvalues. 
We take advantage of these relations in order to ask whether the current and the next generation of short-baseline 
neutrino experiments is capable of definitively ruling out — or, perhaps, ruling "in" — eV-scale right-handed neutrino 
masses for all allowed values of the unknown parameters in U. 

If there are only two right-handed neutrinos (n = 2, fittingly dubbed the minimal model in Q), the situation is quite 
constrained, especially if the neutrino mass hierarchy is inverted. In this case, short baseline neutrino experiments 
should be able to rule out 1714 values less than several eV regardless of the value of 7715 as long as one is sensitive to 
several different oscillation channels. This statement is especially true if one can probe Uf, and disappearance at 
the few percent level, and and I'r appearance at the 10~^ level or better. For 7714 values above 10 eV, the low 
energy seesaw is harder to unambiguously test. A potentially interesting channel to pursue is Vt appearance at the 
10-^ level [3^ or better. 

On the flip side, if light sterile neutrinos do manifest themselves in next generation experiments, it is, in principle, 
possible to determine whether these sterile neutrinos are described by Eq. (jl.ip . Again, in order to do that, one must 
probe several different oscillation channels. In the concrete scenarios spelled out in Eqs. (|III.6IIII.7p . not only did 
we reproduce (close to) the best fit values for the 3+2 fit to current short-baseline data (see [20j), we also "predict" 
that |C/t4,5| is at least as large as |C/^4,5| and |C/e4,5|- This, in turn, would imply that i/^^e ^ Vj- oscillations must be 
observed at a level similar to, or larger than, i/^ -;-> Vf, oscillations. 

For three or more right-handed neutrinos, life is harder. Even if, say, 7774 = 1 eV, 7775 = 10 eV, and TTig = 100 eV, 
it is possible to "hide" the 1/4 state from the active neutrinos at a level such that oscillations mediated by m\ are 



IV. CONCLUDING STATEMENTS AND SUMMARY 
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severely suppressed. Null results can still constrain the parameter space, but the constraints are rather weak even for 
very light right-handed neutrinos, as long as their masses are hierarchical. 

We restricted our discussion to neutrino oscillation experiments. The main reason is that these have the most 
sensitivity to eV-scale sterile neutrinos and, as already summarized above, oscillations allow one to probe different 
combinations of elements of the mixing matrix. For eV sterile neutrinos, however, precision measurements of i3-d eca,y 
spectra and searches for neutrinoless double beta decay (Oj^/3/3) can also provide powerful information [l|, 0, [3^ l34| . 
In both cases, one is sensitive to the Uei elements and the neutrino masses. A detailed estimate of the sensitivity of 
Katrin (tritium /3-decay) can be found in [s^. In the case of life is more interesting. It has been pointed out 

that if all right-handed neutrinos ar e lig ht (masses below a few MeV), Eq. (jl.l[) predicts that the rate for is 
vanishingly small [H, ^ (see also [H, 1351 ]) . A positive signal for 0^/3/3 would definitively exclude the possibility that 
Eq. p.ip is correct and that all right-handed neutrino masses are light. On the other hand, a negative signal for Oiy/3/3 
combined with the discovery of light sterile neutrinos in neutrino oscillation experiments, and the discovery that the 
neutrino mass hierarchy is inverted, would provide a very strong boost for the low-energy seesaw scenario.* 

Finally, the existence of eV scale sterile neutrinos is challenged by cosmol ogy , especially big bang nucleosynthesis 
and the large scale structure of the universe (for recent detailed analyses, see [36, 37]). We have nothing to add to this 
particular discussion. None too pleasing "ways out" have been proposed in in the literature (for many more details, 
see, for example, (38j and references therein). Instead, we choose to view this issue through an optimistic lens: if 
eV scale sterile neutrinos are observed in laboratory experiments, our current understanding of the evolution of the 
universe will also have to change significantly. 

For larger right-handed neutrino masses, the lower bounds on 5,6 discussed here translate into very small upper 
bounds on Uai,5fi- Cosmolog y a nd astrophysics still imply interesting, nontrivial constraints on the low-energy seesaw 
models (see, for example, 0,123) until masses close to 100 keV, and there remains the possibility that one of the 
right-handed neutrinos is the dark matter, as proposed in [sij . as long as all Xaj are small enough for the "dark 
matter" neutrino Vj. Terrestrial experiments, however, lack the sensitivity to unambiguously test the low-energy 
seesaw in case all right-handed neutrino masses are much larger than 1 eV and the lower bounds presented here are 
saturated. Indeed, in order to observe the seesaw neutrinos for such "heavy" right-handed neutrino Majorana mass 
parameters, one needs Nature to pick points in the parameter space such that ^a4,5,6 are much larger than 1 and 
hence ver y fa r from the lower bounds discussed here. This possibility has been considered in the literature (see, for 
example, [i^-lil). 

There are many other consequences of the non-generic properties of U. Before concluding, we highlight the apparent 
double-life of the Majorana phases.^ In the case of two active and two sterile neutrinos, Eq. (jll.24[) expresses the values 
of Ua.b-4:,5 in terms of the neutrino mass eigenvalues, the active mixing angle i?, the complex angle C, and the active 
Majorana phase ip. Measurements of all oscillation channels would reveal, for a generic choice of the parameters, that 
Pab 7^ Pba and the neutrino and antineutrino appearance probabilities would be different, i.e., that CP-invariance 
and that T-invariance are both violated. In turn, we would be able to measure the CP-violating parameters ip and 
^. This would imply that by exploring only lepton number conserving processes one would be able to determine the 
Majorana phase or, equivalently, the Majorana phase is responsible for a lepton-number conserving phenomenon, ^ i.e., 
it determines the magnitude of what is normally referred to as a Dirac phase. While this appears puzzling, it is only 
a reflection of the way we chose to parameterize the mixing matrix and the neutrino masses. Clearly, the definition of 
"input" and "output" parameters is a matter of taste, not physics. Another naive, but much more serious consequence 
one would be tempted to draw from such a measurement is that the neutrinos are Majorana fermions. The logic is 
simple. If the neutrinos were Dirac fermions, the Majorana phase would be unphysical and hence the measurement of 
a Majorana is only possible if the neutrinos are their own antiparticles. The logic is simple, but flawed. One can only 
infer that neutrinos are Majorana fermions because one has made the assumption that the 2-1-2 scenario is correct. 
There is another scenario, with two Dirac active neutrinos and two Dirac sterile neutrinos, that exactly mimics the 
2-1-2 scenario above as far as all lepton-number conserving observables are concerned. It is true that the 2-1-2 seesaw 
would be a much more palatable fit to the data — it would, for example, explain several relations that would appear 
to be accidental in the case of Dirac neutrinos — but the evidence for the Majorana nature of the neutrino would 
remain circumstantial until lepton number violating phenomena were observed. 

We conclude by re-emphasizing the main point of our analysis: the type-I seesaw Lagrangian yields a non-generic 
set of active-sterile oscillation parameters — the masses and mixing angles are, in some sense, entwined. For this 



* Alas, the possibility that the neutrinos are Dirac fermions would, of course, remain. 

t We thank Serguey Petcov for pointing this out at several occasions in the context of leptogenesis and expectations of the minimal 

supersymmetric standard model for charged- lepton flavor violation. See, for example, [43l |. 
-t There is no theorem that prevents this from happening. For a simple concrete example of Majorana phases impacting CP-violating but 

lepton-number conserving phenomena, see [44| . 
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FIG. 5: Allowed values of the magnitude of m^,i as a function of the smallest neutrino mass mughtcst, assuming the known 
oscillation parameters are fixed to Eq. and allowing all phases to vary within their allowed parameter space, for a normal 

(blue, open hatched region) or inverted (green, solid region) neutrino mass hierarchy. The solid lower bounds are proportional 
to the curves presented in Fig. |4l 

reason one is capable to, in principle, test the model by performing enough measurements of neutrino oscillations. 
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Appendix A: Allowed values for the fifi element of the active neutrino mass matrix 

If the right-handed neutrinos are infinitely heavy and integrated out of the theory, one is left with a 3 x 3 "active" 
neutrino mass matrix. Its elements are given by 

mab = ^ {VMNs)ai{VM-Ns)bimi. (A.l) 
2=1,2,3 

Under these circumstances, the ee element of mat, 'm-ee, is the physical quantity probed by Ov^P searches and the 
subject of a lot of theoretical scrutiny (see, for example, [7]). The /i/z element, m^^, plays a similar role for lepton- 
number violating phenomena involving muon-neutrinos, for example, the very rare lepton-number violating kaon decay 
7r=F/i±;i±. For a more detailed discussion see, for example, [s^ . Figure [5] depicts the allowed values of the 
magnitude of m^^ as a function of the smallest neutrino mass mughtestj assuming the known oscillation parameters 
are fixed to Eq. (jIII.l[) and allowing all phases to vary within their allowed parameter space, for both a normal and 
inverted neutrino mass hierarchy. Similar plots for all niab can be seen in [32] . 
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